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TRIGONOMETRIC TIME INTEGRATORS FOR THE 
ZAKHAROV SYSTEM 

SEBASTIAN HERR AND KATHARINA SCHRATZ 


Abstract. The main challenge in the analysis of numerical schemes for the 
Zakharov system originates from the presence of derivatives in the nonlinearity. 
In this paper a new trigonometric time-integration scheme for the Zakharov 
system is constructed and convergence is proved. The time-step restriction is 
independent from a spatial discretization. Numerical experiments confirm the 
findings. 


1. Introduction 

We consider the Zakharov system 

idtE + A E = uE, 
d tt u - A u = A\E\ 2 , 

with initial conditions 

E( 0) = E 0 , it(0) = u 0 , d t u( 0) = ui, 

for given initial data Eq. uq, U\ in appropriate Sobolev spaces. This system is a 
scalar model for Langmuir oscilations in a plasma, see [23, 25]. Here, E : R 1+d —> C 
denotes the (scalar) electric field envelope and u : R 1+d —»■ R the ion density 
fluctuation in spatial dimension d € N. For practical implementation reasons we 
impose periodic boundary conditions, hence both E and u are considered to be 
spatially periodic. 

The Zakharov system has a Hamiltonian structure and conserved quantities. 
More precisely, for strong solutions we have 

4 [ \E(t, x)\ 2 dx = 0 (1.3) 

at Jfd 

and, if u\ has mean zero, 

2 [ \\7E(t,x)\ 2 +u(t,x)\E(t,x)\ 2 + ^\\X7\- 1 d t u(t,x)\ 2 + l-\u(t,x)\ 2 dx = 0, (1.4) 
at Jjd Z Z 

where |V| = \/—A and T d = R d /(2 7rZ) d . The latter is called conservation of energy. 

Several time integrators for solving the Zakharov system numerically have been 
proposed. Due to the outstanding performance of splitting methods for nonlinear 
Schrodinger equations, see the recent papers [9, 10, 20] and references therein, 
splitting methods for the generalized Zakharov system were constructed in [2, 1, 18]. 
In [22] finite differences for the time discretization and a pseudo spectral method 
for the space discretization were used to simulate the Zakharov system numerically. 
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Numerically, the above schemes have been tested extensively. However, due to the 
difficult structure of the system, as explained below in more detail, a convergence 
analysis is missing. 

For the one dimensional Zakharov equations fully-implicit and semi-explicit 
Crank-Nicolson type approximations based on finite difference in time and space 
were derived in [13, 14] and [7, 8], respectively. Numerical experiments [8] indicate 
that the semi-explicit method (which is explicit in n and implicit in E) is preferable 
over the fully implicit method (which is both implicit in n and in E ) due to the 
high computational costs of the latter. However, its convergence only holds under 
the constraint At = Ax, where At and Ax denote the time and space discretization 
parameters. Furthermore, due to the use of the Sobolev embedding theorem the 
convergence results only hold in one dimension. 

The main challenge in the construction and analysis of any numerical scheme 
for the Zakharov system (1.1) originates from the presence of derivatives in the 
nonlinearity: Mild solutions are given by 


E{t) =e ltA E( 0) - i [ e i(t -« )A it(C)F;(^)de, 

Jo , t (1.5) 

u{t) =cos(f|V|)u(0) + s 111 |^M) m /( Q ) + sin (( f _ £)|V|)|V||F;(£)| 2 d£. 


However, it is not obvious how to bound the quadratic term |V||i?| 2 , since “naively” 
estimating the solutions yields 


\\m\\ s < pm+c rMotmom, 

Jo 

ll«(t)lli < \wm\i + 1 |«'( 0)||,_1 + c f iTOIliVid^, 

Jo 


s > d/2, 
1 + 1 > d/2, 


which amounts to a loss of derivatives, see Section 3.1 for a definition of || • || s . 

In order to avoid this, we follow the strategy presented in [21]: We reformulate 
the Zakharov system as a system in (E,dtE,u,dtu). This allows us to construct 
trigonometric time-integration schemes for the Zakharov system (1.1) without im¬ 
posing any spatial-dependent time-step condition or too restrictive regularity as¬ 
sumptions on the initial data (such as analyticity). In particular, their convergence 
also holds in the limit Ax 0. 

For recent developments in trigonometric and exponential integration schemes 
for wave-type equations we refer to [11, 15, 16, 17] and the references therein. For 
local-wellposedness of the Zakharov system in Sobolev spaces of low regularity on 
T d we refer to [6, 24, 19]. Concerning the well-posedness theory on we refer to 
[21, 5, 12, 3, 4] and references therein. 
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2. Trigonometric integrators for the Zakharov system 


To avoid the loss of derivatives we use the method devised in [21]: We reformulate 
the Zakharov system (1.1) as 

r t 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 


id t F + A F = uF + d t u (e{ 0) + F(^)d^j 

dttu — Au = A|£’| 2 , 

(-A + 1 )E = iF-(u- 1) (^(0) + j^ F(Z)dt) , 

where F = dtE (cf. [21]), with initial conditions 

F(0) = i(AE(0) - u(0)E(0)), u(0) = u 0 , d t u{ 0) = u u E{0) = E 0 . 

Let 

X F (t) :=E 0 + [ F{ A)dA. 

Jo 

Then the mild solutions of (2.1) at time t n+ \ = t n + r with to = 0 read 

F(t n + t) =e lTA F(t n ) - i J e^ T -^ A ^uF + u'l F )(t n +0) d£ 

u(t„ + t) =cos(r|V|)u(t n ) + |V| _1 sin(T|V|)u'(t„) 

+ f ivr 1 sin((r - £)|V|)A|£;(t n + £)| 2 d£, 
do 

u'{t n + t) = - |V| sin(r|V|)«,(t n ) + cos(r|V|)u'(t n ) 

+ fcos((T-0|V|)A|£(t n+ 0| 2 d£, 

do 

E(t n + t) =(1 - A) -1 (iF(t n + t) - ( u(t n + t) - 1 )X F (t n + r)). 

In Section 3 we develop a first-order trigonometric integration scheme based on the 
reformulation (2.4) and rigorously carry out its convergence analysis. Furthermore, 
in Section 4 we indicate a generalization to a second-order trigonometric integration 
scheme. 


(2.4) 


3. A FIRST-ORDER SCHEME 

In order to construct a robust first-order scheme we approximate the exact solu¬ 
tions (it, ?/, F, E)(t n +£) appearing in the integrals in (2.4) via Taylor series expan¬ 
sion up to the first-order remainder term. This allows us to integrate e l ^ A , cos(£A) 
and sin(£A) exactly. Furthermore, we use the following approximation for the 
integrals over F : Note that for 0 < £ < r 

/■tn-K 71-1 ftk +1 rt„- )-5 

/ F(A)dA = Y, / F( A)dA + / F(A)dA 

Jo k—0 J^k J tn 

Tl 1 /.t 

= / F(tk + A)dA + / F(t n + A)dA (3-1) 

fc=0 170 ^ 0 

n 

= T^2,F{t k ) + JV,€,n, 

/ c—0 
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where 


^ 1 nT nX /*£ pX 

Fr,£,n ■= (£ - T)F(t n ) + E / / F'(t k + r)drd\ + / / F'(t n + r)drdA. 

fc=n •'° u '° •'0 

(3-2) 


We observe that 


||J>,£,n||s < T||F(t n )|| s +Tt n SUp ||F' (t)\\ 

t€[0,t„+i] 

and it this sense we have, for 0 < £ < r, 

rt n +£ « 

/ F(A)dA«r^F(4). 

•'° fe =0 

Recall the initial conditions (2.2). By setting 
E° = E 0 , u° = u 0 , u° = u l , F° = i(A£° - n°£ 0 ), S° f = E 0 + tF° , (3.3) 
we obtain, for n > 0, the hrst-oder trigonometric time-integration scheme 
1 - e irA 

I „ . I </ Li 1 1 1- I „.'<(■ _ 


^in+l _ girA^ir 


irA 


■ (u n F n + n' 11 ^) . 


w n+1 = cos(r|V|)w n + IVI- 1 sin(T|V|X n + tM" 1 -— C °^| |Vl) A\E n \ 2 , 


u' n+1 = 


= —|V| sin(r|V|K + cos(r|V|)w m + r sm M V D A \ E n\^ 

TV 


(3.4) 


^n+1 _ Qn 


= Sp + rF n+1 , 

E n+1 = (-A + l)" 1 (iF n+1 - (-.u n+1 - 1)S" +1 ) , 

Remark 3.1. Note that for given (E n , F n , u n , u ' n , Sp) we can compute the next 
iteration without saving (E k , F k , u k , u ' k , S F ) for any k < n. 


Remark 3.2. For initial data of sufficiently high Sobolev regularity we will prove 
that the scheme (3.4) is of first-order. Note that one can also use higher order 
quadrature formulas to generate higher order schemes, given additional smoothness 
of the initial data. We give a generalization to a second-order scheme in Section 4. 


3.1. Error analysis. In this section we carry out the error analysis of the trigono¬ 
metric time-integration scheme (3.4). In the following we set for f(x) = Yhkei. d f(k)e lk 
and s£l 

IV|7(*) “ E l*IV(*)e ifc - x , (V>7(*) := |V| °f(x) + /(0) 
feez d 


and define 

11/11.:= ||<V)7IIl W 

For s > d/ 2, we will exploit the fact that H s (T d ) is an algebra, with the standard 
product estimate 

\\fg\\s<c\\f\\ s \\g\\s, 

where c only depends on d and s. Furthermore, we denote by C(X) the space of 
bounded linear operators T : X —> X , and sometimes we write ||T|| S instead of 
||T||£( ffS (T<i)) for the sake of brevity. 

In view of the structure of the Zakharov system 

||(£(t),n(t)V(t))ll W := ||E(i)|| 8+2 + \\u(t)\\ a+1 + ||u'(i)|| 8 
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is the natural norm for our error analysis, the auxiliary function F will be measured 
in || • || s then. 

Theorem 3.3. Fix s > d/2 and 0 < 7 < 1. For any T £ (0, 00), suppose that 
E £ C([0,T];ff s+2+27 (T rf )), u £ C([0, T]; H s+1+21 (T d )) n C 1 ([0, T]; H s+2 "*(T d )) 
is a mild solution of ( 1 . 1 ) with 

m s+2l (T):= sup \\(E(t),u(t),u'(t))\\ [s+2 ^ < 00 . (3.5) 

te[o,T] 

Then, there exists tq > 0 such that for all 0 < t < tq and t n = nr < T the 
trigonometric time-integration scheme (3.4) is convergent of order 7, i.e., 

II (E(t n ) - E n , u(t„) - u n ,u'(t n ) - u ,n )||[ s ] < e Cl c 2 T 7 , 

where ci and c 2 depend only on m s (T) and m s+2l (T), respectively, as well as on 
T, d and s. 

Remark 3.4. Theorem 3.3 implies first-order convergence in the case 7 = 1 . 
Remark 3.5. Note that the Zakharov system (1.1) is locally well-posed in the space 
H s (T d ) x H e ( T d ) x H e -\T d ) 9 (E,u,u'), 

provided that 

0 < s — £ < 1, l/2<£+l/2<2 s, 

0 < s — £ < 1, 1 < ^ + 1 < 2 s, 

0 < s — £ < 1 , d- 1 < t + d/2 < 2 s, 

see [19, 24]. Hence, for 0 < 7 < 1 and 

||-E ! (0)|U+2+27 + ||u(0)|| s+ i + 2 7 + ||m / (0)|| s+ 2 7 < M 1 

there exists a To = To(M 7 ) > 0 such that 

E £ C([0,T o ];ff s+2+27 (T d )), u £ C([0, T 0 ]; 17 s+1+27 (T rf ))nC 1 ([0,T 0 ]; H s+2 ^{T d )), 
which implies that (3.5) holds at least for T = Tq. 

Proof of Theorem 3.3. Let s > d/2. Due to the fact that H s (ff d ) is an algebra it is 
easy to see that the mild solution (E, u) of (1.1) satisfies d t E £ C([0, T\; H s+2 ' 1 (T d )) 
and that (F, E, u) solves (2.4) for F = dtE , see above. In the following, c denotes a 
generic constant which depends on d and s only. We will prove the claim for n + 1 
instead of n. Subtracting the numerical solutions (3.4) from the exact solutions 
(2.4) yields 

F(t n+ 1 ) - F n+1 = e iTA (F(t n ) - F n ) 

1 _ e ir A . 

+ IT \u{t n ){F(t n ) - F n ) + {u(t n ) - u n )F n 

n 

+ ( u'(t n ) - u' n )E{ 0) + u\t n )(Tj2(F(tk) - F k )) 

k =0 
n 

+ (u , (t n )-u ln )(TY / F k ))+L n F , 

k—0 


for d = 1, 

for d = 2. ( 3 . 6 ) 

for d> 3, 


(3.7) 
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and 


<V)(u(t n+1 ) - u n+L ) = cos(r|V|)(V)(u(t n ) - u") 

+ sin(r|V|) j^j-(u'(f n ) - u' n ) 


(3.8) 


as well as 


r|V| 

+ <V)LS, 


u'{t n+1 ) - u' n+ = - sin(T|V|)-^y(V)(u(t„) - u n ) 

+ cos(r|V|)(u'(t rl ) - u' n ) 
sin(r|V|) 


(3.9) 


+ T- 


r|V| 


A (\E(t n )\ 2 — \E n \ 2 ) + L™,, 


and 


(3.10) 


E(t n+1 ) - E n+1 = (—A + l)" 1 (»(*■(*„+!) - F n+1 ) 

n 

- (u(t n+1 ) - u n+1 )(E(0) + rJ2F k+1 )) 

k =0 
n 

+ (1 - u(t n+ i))(r^(F(4+i) - F fc+1 )) + A LI). 

k—O 

The local errors at time t n satisfy 

II^IU = | J (u(t n + $F(t n + 0 - u{t n )F{t n ) 

n 

+ u'(t n +£)I F (t n +£) - u'(t n )(E( 0) + r^F(tfe))^d£ , 

k =0 S 

H<V)iS||. = || jf sin (( T - £)M) A (l^(*r. + Ol 2 - l^n)| 2 ) dc|| s , (3.11) 
" cos((r - 0|V|)A (| E(t n + £)| 2 - \E{t n )\ 2 ) d£ 


\K'h = 


tn+T 


l|A^|| s = (l-u(t n+1 ))( F(X)dX-rY / F(t k+1 )) 

KJo k =o 

By Lemma 3.6 below we have 

max {||L^|| s +||<V)^|| s + ||L^|U+t||AL||U} < C T 1+ ^(l+m s+27 (T)) 4 . (3.12) 

0 <k<n 

Hence, the local errors (3.11) are of order r 1+7 . 

In order to deduce convergence of order 7 globally from (3.12) we need to analyze 
the stability of the integration scheme (3.4). In the following we set 

< = max {||£ fc || s+2 + \\F k \\ s + || U fc || s+1 }. 

11 /C \ 

(i) Error in F: Note that for all s £ R 

l|e iTA |U<l, ||(*rA) _1 (l — e iTA )|| s < 2. 


(3.13) 
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Plugging the stability bound (3.13) into the error recursion (3.7) for F yields that 

| 


\\F(t. n+ i) - F n+1 \\ s < {l + TctnUisitn)) max \\F(tk) — -F fc || s 

0 <k<n 

+ c(m s ( 0) + t n m n s ) (T\\u(t n ) - u"|| fl + T\\u'(t n ) - i/"||«) + \\L n F \\ t 
(ii) Error in ((V)u, u'): We define the operator 

0 _( cos ( r l v l) sin 0lv|)^j-\ 

\-sin(r|V|)|^ cos(t|V|) ) 

Formulas (3.9) and (3.8) imply that 


(3.14) 


(3.15) 


^V)( M (t n+ i)- W n+1 ) N 


u'(t n+ 1) - u' n+1 

+ T 


= O r 


'{v){u(t n )-u n )^ 

u'{t n ) - u’ n 


1—cos(t|V|) (V) \ , 

l,Lwvl) IVI ) a (l £ ( f »)l 2 -| £ "l 2 )+' ■ 

r I v / 


(3.16) 


V T n 

\ L u' 


(3.17) 


Note that the error recursion in E given in (3.10) yields that 
\\E(t n ) - E n \\ s+2 < (1 + ct n m s (t n )) max ||F(4) - F fe || s 

0 <k<n 

+ c(m s ( 0) + t n m™) || u{t n ) - w n || s + || AL” _1 || S , 

which allows us to solve the error recursion in ((V)it, i/) as no loss of derivative 
occurs. More precisely, we have by (3.10) that 

V“(£(i n ) - E n ) = ~ F n ) - ( u(t n ) - u n ) V (t n ) 

n 

+ (1 - u(t n )){rY,{F(t k ) - F k )) + A L^- 1 ) (3.18) 


fc=l 


V| 2 + 
v 


Y v(tn)(V) 1 ((V)(w(t n ) — u n )) + r^ r 


where, for any a = 0,1, 2, || | vp+i IU — 1 and 

IK, JU < (! + ct n m s (t n )) max \\F{t k ) - F fc || s + ||AL^" 1 || i 

0<k<n 

n 

||??(in)|| s = ||-E(0) + r^F J || s < m s ( 0) + t n m™. 
l=i 

Note that for all 0 < r < 1 we have 


sin(r V|) 

< 1, 

1 — cos(t V ) (V) 

r|V| 

s 

r|V| |V| 


< 2 


(3.19) 


(3.20) 


and 


A (\E(t n )\ 2 - \E n \ 2 ) =Re{(E(t n ) + E n )A(E(t n ) - E n ) 

+ (A E(t n ) + A E n )(E{t n ) - E n ) 

+ 2V(£(f„) + £") • V(£7(t n ) - E n )}. 


(3.21) 
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In particular, 

|| A (\E(t n )\ 2 - \E n \ 2 ) |U < Ac\\E{t n ) + E n \\ s+2 \\E(t n ) - E n \\ s+2 
< 4c(m s (i„) + to") \\E(t n ) - E n || s+2 . 
Plugging (3.20), (3.21), (3.18) and (3.19) into (3.16) we obtain that 
/(V)(u(i n+1 ) - u" +1 )^ =(n ^ ((V)(u(t n ) - u") 

V u'(t n+ 1 ) - u' n+l 


— u 


where 


— (Or + P T ,n) I 

\ u'(t n ) 

/ 1—cos(r]V|) <V> n \ 

P r,k=T\ ^|v|)' V| n )M Pk+Pk+Pk}, 

V r|V| U / 


+ E-T,ni (3.22) 


with the operators 


Pk = ( E(t k ) + E k ) A(|V| 2 + 1)" Vi fc )(V>-\ 


(3.23) 


(3.24) 


p 2 fc = (A£(t fc ) + A£i fc )(|V| 2 + l)- 1 r 7 (t fc )(V>- 1 , 
p 3 k = 2(\7E(t k ) + \7E k ) • V(|V| 2 + ^“^(tfcKV)- 1 . 

The remainder satisfies 

llftr.nll. < c|| E(t n ) + E n \\ s+2 \\r 2 T J s + ||<V)LS|| 8 + ||T"/1|s 

< CT(m s (t n ) + m")((l +t n m s (t n )) max \\F(t k ) - -F fe || s + ||AL^ _1 || S ) 

0 <k<n 


+ \\(V)LZ\\ S + \\L”,\\ S . 


(3.25) 


Note that (3.19) implies that for all / G H a 1 (T d ) and j = 1, 2, 3 we have 
max \\p J k f\\ s < c max \\E(t k ) + E k \\ s+2 \\r)(t k )\\ s \\f\\ s - 1 

0<k<n 0<k<n /o o^\ 

< c(m s (t n ) + m")(m s ( 0 ) + t n m")||/|| s _i. 

Thus, plugging (3.20) and (3.26) into (3.15) we obtain that 

max sup ||T _ 1 P r , fc || £((H « (T d )) 2 ) < c(l + t n )(m s (t n ) + m n s f =: q n , ( 3 . 27 ) 
U</c<n o<r<ro v 

The bound (3.27) and Lemma 3.7 below yield the essential stability bound 


P (O t + P T ,k) | 

k—ko 


C{{H «(T<*))2) - 


< e 1 


*(l+9~) 


(3.28) 


for any 1 < fco < n. Thus, solving the error recursion in (3.22) we obtain by the 
stability bound (3.28) and the bound on lZ T , k in (3.25) that 

\\u(t n+1 ) - u" +1 || s+ i + ||u'(t n+ i) - u'" +1 || s 

n 

^ 2 n max ||7^. r ^|| s max II I I (O r H- -Pr,fc)|| r((H s (T d ^'\ 2 ^ 

0 <k<n l<ko<n AJ - M (il )) ) 

- ~ “ k—ko 

(ctn(m s (tn) + to "){(1 + t n m s (t n )) max \\F(t k ) - -F fc || s 


(3.29) 


< 


+ max 

0 <fc<n 


ax HAillU} +2n max (||(V}T*|| S + ||Lj;,|| a ))e t »( 1 +«-). 

c<n 0<fc<n / 
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The error bounds (3.14), (3.17) and (3.29) together with the bound on the local 
errors in (3.12) yield 


||F(t n+1 ) - F n+1 \\ s < (1 + ^(m?)) max || F(t k ) - F k \\ s + A 2 (3.30) 

0 <k<n 

and 


\\E(t n+1 ) - E n+1 \\ s+2 < Ai(m” +1 ) max ||F(f fc ) 

0<fc<n+l 


F fc || s + t 1 A 


2 , 


\\u(t n+1 ) - M n+1 || s+ l + \\u'(t n+1 ) - U ,n+1 \\s 
< Ai(m”) max \\F(t k ) - F k \\ s + t 1 A 2 , 

0 <k<n 


(3.31) 


(3.32) 


where A\ = Ai(-) is a continuous and monotonically increasing function which also 
depends on m s (T), A 2 is a constant which depends on m s + 2j (T), and both A\ and 
A 2 depend on T, d and s. From (3.30) we obtain 

n 

max \\F(t k ) - F fc || s < A 2 V (l + t^M) V+^ < Te TA ^A 2 rF (3.33) 

0<fc<n+l v ' 

- 3=0 

Then, (3.32) implies 

max {||u(t fc ) - W fc || s+ i + \\u'(t k ) - v! k 11 s } < ( Te TA ^ + 1)A 2 tT (3.34) 

0<fc<n+l 

Similarly, (3.31) implies 

max || E(t h ) - E% +2 < {A^m^Te™^ 1 ^ + 1 )A 2 t\ (3.35) 

0<fc<n+l 

Now, the assertion follows by a continuity argument: We obtain that 

m n s +1 < 2( J 4 1 (m? +1 )Te TAl « +1) + 1 )A 2 t^ + m s (T) 

The quantity m™ +1 depends continuously on r and tends to zero as r —> 0. We 
conclude that m” +1 < 2 m s (T) as long as 

0 < r < m s (T)^{2(A 1 (2m s (T))Te TAl ^ m ‘ ( - T '>'> + 1 )A 2 )~$ =: r 0 . 

The claimed estimate (for n+ 1) follows with the constants c 2 = 4TAi(2m s (T))7l 2 
and Ci = TAi(2m s (T)). □ 

Lemma 3.6 (Local error). Let s > d/2. For 0 < 7 < 1 the local errors defined in 
(3.11) satisfy 

max {||L$.|| 8 + ||<V)l£|| a + ||L^|| S +r||AL||| s } < CT 1+ ^t n (l + ro 5+27 (T)) 4 , 

0<k<n 

where m s + 2 "/(T) is defined in (3.5) and c depends on d and s. 

Proof. In the following fix 0 < 7 < 1 and let c denote a constant depending on s 
and d only. 
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The mild formulations (1.5) and (2.4) yield 
E(t n + £) — E(t n ) 

e i?A _ 1 r t 


(-£A)7 

F(t n + £) — F(t n ) 

e i«A _ x 


(-£A yE(t n )-i [ e^-^ A (uE){t n + A)dA, 
Jo 


(-£A)7 


(-£A) 7 .F(t„) - i [ e^-^iuF + u'E){t n + A)dA, 
Jo 


and 


|V|(u(i n +0 - u{t n )) 

s C ° Ki'vlr" + 

+ / S sin((£ — A)|V|)A|£ , (t ra + £)| 2 , 
Jo 

u'{t n +£) - u'(t n ) 

+ cos((£ - A)|V|)A|£7(t„ + Ol 2 - 
Jo 


Note that 


sin(^|V|) 




< i, 


1 — cos(£|V|) 


< 2 , 


1 - e^ A 


(~W 


< 2 . 


(£|V|) 7 

Plugging (3.38) into (3.37) yields for 0 < £ < 1 

|| u(tn + £) u(tn) || s+1 + ||u'(tn + 0 — u(t n )\\ s 

< c£ 7 (||u(i„)||s+ 1+7 + |K(*n)||a+'r) + c£m s (T) 2 

< c£ 7 (l + to s+7 (T)) 2 . 

We have 

ll( — A) 7 F(t)|| s = ||(—A) 7 .E'(f)|| s 

<||(-A) 7 + 1 £(t)|| s + c||(-A)Mt)|| s ||(-A) 7 J E;(t)|| s . 


Hence, plugging (3.38) and (3.40) into (3.36) yields for 0 < £ < 1 
II E(t n + £) — E(t n )\\ s+2 + || F(t n + £) — i ? (tn)||s 

< cC\\E(t n )\\ s+2+ 2 j + c C0 + m s (T)) 2 + c£(l + m s (T)) 3 

< c£ 7 ( 1 +m s+ 2 7 (T)) 3 . 

(i) Local errors (V)L™ and L By the definition of (V)L™ and L ”, 
we obtain with the aid of (3.41) that 

||<V)X”|| s + ||L”,|U<ct sup ||A(|J57(£„ -hOl 2 - l^(in)| 2 )|| s 

0 <f<T 

< crm s (T) sup \\E(t n + £) - E(t n )\\ s+2 < cr 1+7 (l + m s+27 (T)) 4 . 
o <£<r 


(3.36) 


(3.37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 
in (3.11) 


(3.42) 
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(in) Local errors Lp and A Lp: By the definition of A Lp and L ^ in (3.11) we 
have 

rt n +T n 

\\AL n E \\ s < cm s (T)\\ / F(A)dA-r^F(t fc+ i)|| s , 

"'° fc=o 


j f\\s < C J m s(T)(\\F(t n + (;) ^(tn)||s + ll-Tv^ralls^ (3 43) 

+ (1 + m s (T)) 2 (||u(f„ + 0 - u{t n ) || s 
+ (1 + t n )\\u'(t n + £,) — u , (t n )||s^d^, 

cf. (3.1) for the definition of P T ,^, n - Note that by (3.36) we obtain for 0 < £ < r 
that 

rt n +£ n 

J r T,t,n= / F(A)dA-T^F(ifc) 

"'° fc=0 

n—1 « T />£ 

= / (F{tk + A) — i 7 ’(t/ e ))dA + / -F(£ n + A)dA — rF(t n ) 

n </o -'O 


fc =0 
n—1 


r iAA _ i 

E/ ^TT7(-^ A ) 7 ^(4)dA + ^„, 


fc =0 ’ 


(-AA)i 


where 


"-1 /.t 


nT nA 

HEJU < II E / / e i(A -« )A (uF + u'E)(t n + £)d£dA|| s + rcm s (T) 2 

k= o"' 0 •'° 

< Tct n (l + m s (T)) 3 . 

Hence, (3.38) together with (3.40) implies that for 0 < £ < r 

||Jv,i,n||a < cr 7 t n m s+ 2 7 (T) + cr 7 i n (l + m s (T)) 3 . 

Plugging (3.39), (3.41) and (3.45) into (3.43) we obtain that 

t\\AL^\\ s + \\L n F \\ s < cr 1+7 t n ( 1 + m s+27 (T)) 4 . 
Collecting the results in (3.42) and (3.46) yields the assertion. 

Recall the definition of the operator O t from (3.15). 


(3.44) 


(3.45) 

(3.46) 

□ 


Lemma 3.7 (Stability lemma). Let s€R. For 0 < r < tq, \ < ko < k < n let 
P T ,k 6 £((.fP(T d ) 2 ) such that 

q:= max sup ||r _1 P rife || £ ( (ffS ( T d))2) < oo. 

0<T<To 

Then, for all ( f,g) € (H s (T d )) 2 


II (°r + P r-k) 

k=ko 



< e nr(l+q) 


W 

s 

W 


Proof. Let 


V = 


J_ (l i 

V2V 1 


, Z T 


0 -^«V)-|V|) A 

, 2A ^«V)-|V|) 0 , 
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We have 

O t = V _1 diag (e iTlvl ,e~ iTlvl )V + Z T . 

Note that the action of Z T is nothing but multiplication by r of the zero mode of 
the second component. For Q T ^ = V(Z T + P Tjfc )I / _1 we obtain 


n n 

n (Or + P T , k ) = V~ 1 { [] (diag(e^l v l, e-^' v l) + }w. 


k=ko 


k=ko 


Hence, 


(Or + Pr,k) 


< n (l + IIQr,fe||z;((ff*(T< J )) 2 ) 

k—ko / k—ko 

Due to 

IIQ-r,fc||£((i/o(T d )) 2 ) < II Z T + -P-r,fc||,C((ff s (T' i )) 2 ) < T + Tq 
we conclude that 


(3.47) 


n 

(l + |IQr,fe||£((iI=(T< J )) 2 )) < (l + r (l + <?)) <e nT ^ 1+q \ 


k—ko 


and the claim follows from (3.47). 


□ 


3.2. Error analysis for strong solutions and in the energy space. 

Remark 3.8. As lower order Sobolev norms are controlled by higher order Sobolev 
norms Theorem 3.3 also yields a convergence result for strong solutions (i.e., in 
H 2 (T d ) x i7 1 (T d ) x L 2 (Y d )) as well as in the energy space (i.e., in H 1 ( T d ) x L 2 (T d ) x 
H~ 1 { T d )). More precisely, assume that for some 7 > 0 the regularity assumptions 
(3.5) hold with s = d/2 + e for any e > 0. Then there exists a To > 0 such that for 
all 0 < r < To and t n = nr < T the following convergence bounds hold: 

\\(E(t n ) - E n ,u(t n ) - u n ,u'(t n ) - u' n )\\[ r ] <ct 7 , r = —1,0. (3.48) 

However, the regularity assumptions on the data are quite strong. 

In the following we will show that in dimensions d < 3 the regularity assumptions 
(3.5) with s = max(l,<i/2 + e) actually imply first-order convergence , i.e., (3.48) 
holds with 7 = 1 . Here, we apply asymmetric product estimates and in order to 
control the error of F and u' in L 2 (T d ) and H~ 1 (T d ), respectively, we need a priori 
bounds on the numerical solutions in higher order Sobolev spaces, cf. [ 11 , 20 ], 

We will carry out the error analysis in detail only for the energy space as the 
result for strong solutions follows along the same lines. Furthermore, for the sake 
of clarity of the exposition, we restrict ourselves to dimensions d < 3, where the 
following product estimates are crucial for our analysis: For si + S 2 > 0 and 1 < 
d < 3 we have 

\\fg\\s < c\\f\\ Sl \\g\\ S2 for all s < si + s 2 - f with s 1; s 2 and - s ^ f 
\\fg\\s <c\\f\\ Sl \\g\\ S2 for all s < si + s 2 - | with si, s 2 or - s = | 
such that in particular we obtain for 1 < d < 3 and e > 0 that 

ll/sll-l < c||/||—l||s||max(d/2+ e ,l)- 


(3.50) 
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Theorem 3.9. Fix 1 < d < 3 and 7 > 0. For any T £ (0,oo) and any e > 0, 
suppose that for S := max(d /2 + e, 1 ) 

E £ C{[ 0, T]; H 2+8+21 (T d )), u £ C([0, T\,H 1+s+2l ( T d )) n C\[0, T); H 5+2ry (T d )) 
is a mild solution of ( 1 . 1 ) with 

mi+ 2 7 (T) := sup \\(E(t), u(t), u'(t))\\ [ 5+27 ] < 00 . (3.51) 

*£[0,T] 

Then, there exists tq > 0 such that for all 0 < t < tq and t n = nr < T the 
trigonometric time-integration scheme (3.4) is first-order convergent in the energy 
space, i.e., 

| \(E(t n ) - E n ,u(tn ) - u n ,u'(tn) - U ,n )||[_i] < CT, 
where c depends only on ms(T), T and d. 

Proof. In this proof we proceed similarly to the proof of Theorem 3.3. However, 
we need to be more careful when estimating the nonlinear terms. In the following 
fix 1 < d < 3, e ,7 > 0 and set S = max(d/2 + e, 1). First note that the regularity 
assumptions (3.51) together with Theorem 3.3 (choosing s = S) imply that there 
exists a To > 0 such that for all 0 < r < To and t n = nr < T we have 

:= max {\\E k \\ 2+g + \\F k \\ s + ||u fe ||i +5 + ||-u ,fc |U} < 2 m s {T) < 00 . (3.52) 

0 <k<n 

In the following we assume that t < tq such that (3.52) holds. Furthermore, 
we denote by c a constant depending only on ms(T), T , d and prove the claim for 
n + 1 instead of n. 

The regularity assumptions (3.51) imply that the local errors defined in (3.11) 
satisfy 

max {\\L%U + ||(V )^||_ 1 + H^'IU + tWAL^} < ct 2 , (3.53) 

0 <fc<n 

see Lemma 3.11 below. In order to deduce first-order convergence globally from 
(3.53) we need to analyze the stability of the integration scheme (3.4) in the energy 
space. 

(i) Error in F in H~ l : The error recursion in (3.7) together with the stability 
bound (3.13), the bilinear estimate (3.50) and the local error bound (3.53) yields 
that 

\\F(t n +i) - F n+1 ||_i < (1 + Tct n m s (t n )) max || F(t k ) - F fe ||_i 

0 <k<n 

+ c(m s {0) +t n max \\F k \\ s )\\u'(t n ) - u ,n ||_i + c||F rl ||, 5 ||'u(t n ) - u n ||i + ct 2 . 

0 <k<n 

Furthermore, the a priori boundedness of the numerical solutions (3.52) implies 
that maxo<fc<n ||^ fc ||(5 < 2 ms(t n ) < 00 . Hence, we obtain that 

||F(t„ +1 )-F”+ 1 ||_ 1 < (1 + tc) max \\F(t k ) - F k \U 

°<k<n ( 3 . 5 4 ) 

+ Tc(\\u'(t n ) - U n || _1 + || U(t n ) - U"||_i) + CT 2 . 

(ii) Error in E in FL 1 : Similarly we obtain by the error recursion (3.10) together 
with the bilinear estimate (3.50), the bound on the numerical solutions (3.52) and 
the local error bound (3.53) that 

\\E(t n+1 ) - F n+1 ||i < c max ||F( 4 ) - F fc ||_i + c\\u(t n ) - u n ||_i + ct. (3 55 ) 

U \ /C \ 77, 
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(Hi) Error in (u,u') measured in L 2 x H~ x : The a priori boundedness of the 
numerical solutions (3.52) together with the bilinear estimate (3.50) implies that 
for 0 < a < 2 we have 


||(|V| 2 - Q ^(t„) + \X7\ 2 ~ a E n )\X7\ a (E(t n ) - E n )||_! 

< c\\E(t n ) + E n \\ 2+s \\E(t n ) - E n ||r < 2cm s {t n )\\E(t n ) - E n \\ 1 . 


Thus, similarly to (3.22) we obtain that 


Av)(«(t„ +1 )- U " +1 )\ 
\ u'{t n+ 1 ) - u’ n+1 J 


(O t + P T ,n) 


f(V)(u(t n ) - u")\ 
\ u’(t n ) - u ,n J 


E-T,r, 


where 


(3.56) 


ll^r,n||-i < cr( max ||F(t fc )-i ;,fc ||-i + ||AL3|||_ 1 ) + |KV)LS||_i + ||LS,||_i (3.57) 

0 <k<n 

and O r , P Ttn are defined in (3.15) and (3.23), respectively. The bilinear estimate 
(3.50) together with the a priori boundedness of the numerical solutions (3.52) and 
the definition of ??(tfe) in (3.19) furthermore implies that 

||(A E(t k ) + AE k )rj(tk)\\-i < c\\E(t k ) + E k \\ 1 \\ V (t k )\\ s < 2c(l + t n )m 2 s {T). 

Thus, by the definition of p* k in (3.24) we have for all / £ H~ 1 ( T d ) that 

max max [|p^/||_i < c||/||_i. 

j —1,2,3 0<fc<n 

Together with (3.20) (which holds for all s £ R) this yields by the definition of P T:k 
in (3.23) that 

max sup Wt - 1 P Tt k\\c((H- 1 n d )) 2 ) < c - 

0<k<n 0<t<To 

Hence, solving the error recursion in (3.56) we obtain with the aid of the stability 
Lemma 3.7, the bound on 7 Z TiTl given in (3.57) together with the local error bound 
(3.53) that 

IKWi) - u n+1 \\ 0 + |K(i n+ i) - u m+1 ||_i 

<c max \\F(tk) - F k \\- 1 + CT. (3.58) 

0 <k<n 

Collecting the results in (3.54), (3.55) and (3.58) yields the assertion. □ 

Remark 3.10. Note that in the limit t —> 0 Theorem 3.9 (together with Remark 
3.5) implies first order-convergence in the energy space if for some e > 0 

IK-E’(O), u(0), u , (0))||[i +e ] < oo ford = 1 , 2 , 

IK-E’(O), u(0), u , ( 0 ))||[ 3 / 2 +e] < oo ford = 3. 

Lemma 3.11 (Local error in the energy space). Let 1 < d < 3. Then the local 
errors defined in (3.11) satisfy for any e > 0 and S > max(d/2 + e, 1) 

max {||L* ||_r + ||(V)L*|U + H^ll-t + r||AL* ||_r} < cr 2 (l + m s {T))\ 

0 <k<n 

where m$(T) is defined in (3.51) and c depends on T and d. 
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Proof. The strategy of proof is similar to the one of Lemma 3.6. In the following 
fix e > 0 and set S = max(d/2 + e, 1). Let c denote a constant depending only on 
d . The local error representation (3.11) together with the bilinear estimate (3.50) 
implies that 

||^||-i < cf m s {T)(\\F(t n +0- F(t n ) H-! + ||■F T ,*, n ||_ 1 ) 

+ (1 + ms(T)) 2 ^||u(f„ + £) - u{t n ) ||_i 

+ (1 + t n )\\u'(t n + £) — it'(t n )||_i^d£, ^ 

ll(V)I'"||- 1 + ||^||-i < cms(T) f || E(t n + 0 - E(t n ) ||id^, 

Jo 

||AL’|||_ 1 < C m i (r)|| / F(A)dA-r^F( 4 + i)||_i. 

, ' 0 k=0 

Choosing 7 = 1 in (3.37) we obtain with the aid of (3.38) (which holds for all s £ R) 
and the bilinear estimate (3.50) that 

II u(t n + £) - u(tn)||o + | \u'(t n + £) - u'(t n )||_1 

< c^(||u(t„)||i + ||u , (t„)||o) + cf,m s (T) 2 (3.60) 

< c£(l + m$(T)) 2 . 

Note that 

||(-A)F(t )||_ 1 < \\F(t)\\i = P'Wlli < PWlls + cm s (T) 2 . (3.61) 

Choosing 7 = 1 in (3.36) we thus obtain by (3.38), (3.50) and (3.61) that 
|| E{tn + 0 - E{tn) 111 + || F(t n +Q- F(t n ) ||_1 

< c£||£(i n )|| 3 + c£(l + m s (T)) 3 (3.62) 

< c£(l + ms(T)) 3 . 

Similarly we can show that for 0 < ( < r we have (see also (3.45)) 
rtn+i 71 

II / F{\)&\-T^2,F(t k )\\ s <cTt n (l + ms(T)) 3 . (3.63) 

170 k—0 

Plugging (3.60), (3.62) and (3.63) into (3.59) yields the assertion. □ 

For strong solutions we obtain the following convergence result: 

Theorem 3.12. Fix 1 < d < 3. For any T £ (0,oo), suppose that 

E £ C{[ 0, T]; H 4 (T d )), u £ C([0, T]; H 3 (T d )) n C 1 ([0, T]; H 2 (T d )) 
is a mild solution of ( 1 . 1 ) with 

m 2 (T) := sup \\(E(t), u(t), u'(f))|| [ 2 ] < 00 . (3.64) 

te[o,T] 

Then, there exists tq > 0 such that for all 0 < t < To and t n = nr < T the scheme 
(3.4) is first-order convergent in the sense that 

\\(E(tn) - E n ,u(tn) - u n ,u(tn) - m'")||[ 0 ] < CT, 
where c depends only on m 2 (T), T and d. 
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Proof. Fix 0 < e < 1. Note that the regularity assumptions (3.64) imply that 
there exists a To > 0 such that for all 0 < r < To and t n = nr < T the numerical 
solutions satisfy for S = d/2 + e the a priori bound 


m n 5 := max {\\E k \\ s+2 + ||F fc || 5 + \\u k \\ s+1 + ||n ,fc |U} < 2m 2 (T) < oo. 

0 <k<n 


This follows from choosing s = d/2 + e and 7 = 1 — s/2 in Theorem 3.3, whereupon 
in particular 7 = 1 — d /4 — e/2 > 0 as d < 3. Together with the L 2 (T d ) estimate 


||/$||o<c||/||o|M|i, 


(3.65) 


the proof can be completed along the lines of the proof of Theorem 3.9. 


□ 


4. Second-order Scheme 

In this section we derive a second-order trigonometric integration scheme for the 
Zakharov system (1.1) based on the mild solutions (2.4). In order to achieve this we 
use a second-order exponential integrator in the approximation of F. Furthermore, 
we approximate the integrals in (u,u') with a trapezoidal rule, i.e., we use that 



Uf(r)+fm + Os(r 3 sup wnob), 

Z 0<£<T 


(4.1) 


where for notational simplicity we here use the notation O s (z) which denotes a 
remainder term depending on z > 0 when measured in H s , i.e., 


f = g + O s {z) if \\f-g\\ s <cz, 


for some constant c > 0 . 

Using the second-order Taylor series expansion 


(uF + u'X F ){t n + 0 = (uF + u'l F )(t n ) + £(uF + u'l F )'(t n ) 
+ O s (t 2 sup \\{uF + u'Z F )"{t n +£)|| s ) 

V 0<f<T > 


in the integral in F as well as the trapezoidal rule (4.1) for the approximation of 
the integrals in u and v! yields the following approximation to the solutions (2.4) 
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of the Zakharov system: For sufficiently smooth solutions we have 

F(t n + r) =e lTA F{t n ) - i f e^-^d^uF + u'l F )(t n ) 

Jo 

-i [ e^-^d {(u'F + uF 1 + u"l F +u'l' F )(t n ) 

Jo 

+ O s (t 3 sup ||(uF + u , I P ’)"(i n +£)|| s ') 

V 0 <£<T ' 

u{t n + t) =cos(r|V|)u(i n ) + |V | _1 sin(r|V|)u , (t Tl ) + ^ A\E(t n )\ 2 

/ o \ (4.2) 

+ o s (t 3 sup |||£;|"(t„ + 0IU+i) 

V 0 <{<T ' 

u’(t n + t) = - |V| sin(r|V|)u(t„) + cos(r|V|)u'(i n ) 

+ ^ (A| E(t n + t)| 2 + cos(T|V|)A|^(t ri )| 2 ) 

+ O s (t 3 sup +OIU+ 2 ) 

V 0 <f<r ' 

E(t n + t) =(1 - A) -1 (iF(t n + r) - (u(t n + r) - 1 )l F {t n + t)) . 

In order to derive a robust scheme we integrate the terms involving e ^ A £ 5 with 
<5 = 0,1 exactly, i.e., we will use that 


(4.3) 


= __L (1 _ e *' rA ) /, 

= __L {r + J_ ( !- e —))/. 

Next we need to derive a suitable approximation to If defined in (2.3). We have 

/>£ n n—1 « r 

I F (f n ) = E 0 + F{ A)dA = £7 0 + V / ^(tfc + A)dA. (4.4) 

do fc=0 do 

Note that by (2.4) and (4.3) we have 

pT pT pT p A 

/ F(t k + X)dX = e lXA F(t k )dX-i / e l(A - f)A (wF + w'I F )(4+^dedA 
Jo Jo Jo Jo 

pT pT p A 

= / e iAA F(4)dA - i / / e i(A - f)A (wF + w'I F )(4)d£dA 

Jo Jo Jo 

+ (9s (t 3 sup ||(uF + u'Xp)'(tfc+ £)IU) 

0 <4<r 

= 1 (e- A -l)F(4) 

k 

+ ^ ( T - -^(e irA _ X )) ( U (4)F(4) + u'(t k )(E 0 + F(tj)^j 

3 =0 

+ (9 S (r 3 sup ||(uF + u':Z:f)'(4 + OIU))- 
o<c<r 


(4.5) 
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Plugging the above expansion into (4.4) yields that 

n— 1 

l F (t n ) = Eo- F(t k ) 

k—0 

n— 1 k 

+ tV 2 {t A) ^ (u(tk)F(tk) + u'(t k )(E 0 + Ty^F(tj))) 

k—0 j =0 

+ O s (r 2 t n sup ||(u.F-|-u'Xp’) , (£)|| s ), 


where 


i _ f> irA 

Pi(rA):= . , T> 2 (tA) := A _1 (l +X> 1 (rA)). (4.7) 

irA 

Using the differential equations (2.1) as well as the definition of I F in (2.3) we 
furthermore obtain that 


(u'F + uF' + u"I F + u'l'p) (t n ) 

= ( u'F + u(iAF — iuF — iu'Tp) + I F A(u + \E\ 2 ) + u'F ) (t n ). 


Plugging the relations (4.3), (4.6) and (4.8) into (4.2) yields a second-order trigono¬ 
metric time-integration scheme by setting 

E° = Eq, u° = uo, u'° = u±, 

F° = i(AE° - u°E °), S° F = tF°, 1 f := E 0 


and for n > 0 

pn+l = e irApn + ^ ( T A) (u n F n + U ’ n l £) 

+ tV 2 (tA) (2u ,n F n + iu n (AF n - u n F n - u’ n l JJ) + X£A(u” + |U n | 2 )), 
u n+1 = cos(t| V|)M n + M- 1 sin(r|V|)M m + ^Vl” 1 sin(r|V|)A|^ n | 2 , 

n 

l n F +1 =Eo- D 1 (tA)S f + tV 2 (tA) J 2 (u k F k + u ,k (E 0 + S k F )), 

k—0 

E n+1 = (-A + I)" 1 ( iF n+1 - (u n+1 - 1 )I” +1 ) , 

u' n+1 = — |V| sin(r|V|)u” + cos(t| V|)'u ,n + ^ (A|E" +1 | 2 + cos(T|V|)A|U n | 2 ) , 
S" +1 = + rF n+1 . 

(4.10) 


Remark 4.1 (Second-order convergence). For sufficiently smooth solutions the trigono¬ 
metric integration scheme (4.10) is second-order convergent without imposing any 
spatial-dependent time-step restriction, i.e., also in the limit A* —> 0. More pre¬ 
cisely, Theorem 3.3 holds for (4.10) with 7 = 2 . The ideas in the error analysis 
are thereby similar to the ones used in Section 3.1. The only additional important 
estimate is that 

||T> 2 (rA)((A/) ff )|| s < c||/y|ff|| s+2 

for some constant c > 0. We omit the details of the proof and refer to [11] for the 
analysis of second-order trigonometric integrators for semilinear wave equations and 
to [17] for the analysis of higher-order exponential integrators. 
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Remark 4.2. Note that for given ( E n , F n , u n , u' n , Sp,I F ) we can compute the next 
iteration in (4.10) without saving (E k , F k , u k , u' k , S F ,Ip) for any k < n by setting 

I n F +1 := l n F - tT>i(t A)F n + tV 2 (t A)(u n F n + u' n {E 0 + S£)), 1° F = E 0 . 

(4.11) 


5. Numerical experiments 

In this section we numerically confirm the first-, respectively, second-order con¬ 
vergence rate of the trigonometric time-integration schemes (3.4) and (4.10) towards 
the exact solutions of the Zakharov system (1.1). Furthermore, we numerically test 
the geometric properties of the trigonometric integration schemes, i.e., the conser¬ 
vation of the L 2 norm of E(t) (see (1-3)), the conservation of the energy (see (1.4)) 
as well as the shape preservation of solitary waves over “long times”. 

Remark 5.1. Note that in the derived convergence bounds, the error constants de¬ 
pend on T, which is natural in subcritical regimes. Nevertheless, the numerical find¬ 
ings suggest that for a sufficiently small CFL number the geometric quantities are 
preserved on “long” time intervals. Thereby, we define the value CFL := t(Ax ) -2 
with t and Ax denoting the time- and spatial-step size, respectively. 

Remark 5.2. In the numerical experiments we use a standard Fourier pseudo- 
spectral method for the space discretization. For sufficiently smooth solutions the 
fully discrete error then behaves like r + K~ r for the first-order scheme and like 
t 2 +K~ r for the second-order scheme for some r > 0 depending on the smoothness 
of the solutions. For a fully discrete analysis of exponential-type time integrators 
coupled to a spectral approximation in space for Schrodinger, respectively, semilin- 
ear wave equations we refer to [9] and [11], respectively. 

Example 5.3. We consider the Zakharov system (1.1) set on the one-dimensional 
torus T with initial values 

.E^O, x) = (2 — cos(x) sin(2x )) _1 sin(2x) cos(4x) + i sin(2x) cos(x), 

u( 0 ,x) = (2 — sin( 2 x ) 2 ) -1 sin(x) cos( 2 x), d t u( 0 ,x) = (2 — cos( 2 x ) 2 ) -1 sin(x) 

(5.1) 

normalized in Ft 2 , Ft 1 and L 2 , respectively. In order to test the convergence rate 
of the first-, respectively, second-order trigonometric time-integration scheme (3.4) 
and (4.10) we take the numerical method presented in [1] as a reference solution. 
For the latter we choose a very small time-step size to ensure to be sufficiently 
close to the exact solutions. For the space discretization we choose the largest 
Fourier mode I\ = 2 10 (i.e., the spatial mesh size Ax = 0.0061) and integrate up to 
T = 1. The error of (E,u) measured in the corresponding discrete H 2 x FI 1 norm 
is illustrated in Figure 1. 

Example 5.4 (Solitary waves). Exact solutions of the Zakharov system (1.1) are 
explicitly given by so-called solitary wave solutions, which for the Zakharov system 
set on R are described by 


E(t, x) = V2B 2 (1 - C 2 )sech(B(x - Ct)) exp (i {C/ 2x - (C 2 /4 - B 2 ) t ))), 
u(t, x) = — 2B 2 sech 2 ( B(x — Ct)) , (5.2) 

dtu(t, x) = —4E? 3 C'sinh (B(x — Ct)) cosh -3 ( B(x — Ct)) 
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Figure 1. Orderplot (double logarithmic). Convergence rates of 
the first-order scheme (3.4) (blue, circle) and the second-order 
scheme (4.10) (red, star). Left picture: Error in E measured in 
H 2 . Right picture: Error in u measured in H 1 . The slope of the 
dashed-dotted and dashed line is one and two, respectively. 


with B,C £ I. For the numerical simulations we choose “a large torus” (more 
precisely x G [— 107r, IOt]). In Figure 2 we simulate the soliton solution (5.2) with 
the trigonometric integration schemes (3.4) and (4.10) up to T = 100. We carry 
out the simulations for two different CFL numbers. Furthermore, we set B = 0.5 
and C = 0.15. 


t = 50 t = 50 




Figure 2. Simulation of solitary wave \E n \ at different times t 
with first-order scheme (3.4) (red) and second-order scheme (4.10) 
(yellow). Initial profile: blue. Left picture: CFL= 3.2. Right 
picture: CFL= 32. 


Example 5.5 (Energy conservation). In this example we numerically test the L 2 
conservation (1.3) and the energy conservation (1.4) of the first-and second-order 
trigonometric time-integration scheme (3.4) and (4.10), respectively. The numerical 
findings are illustrated in Figure 3 (left picture: first-order scheme, right picture: 
second-order scheme). In both simulations we choose CFLft: 5. For a too large 
CFL number additional numerical findings suggest that the energy is no longer 
conserved. 








TRIGONOMETRIC TIME INTEGRATORS FOR THE ZAKHAROV SYSTEM 
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Figure 3. Simulation of the deviation of the numerical energy 
H(E n , u n , u' n )—H(E °, u°, u'°) and the L 2 norm ||^ n || L 2 — ||_E°|| L 2 . 
Left picture: First-order scheme (3.4). Right picture: Second-order 
scheme (4.10). 
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